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Abstract 

Many articles have addressed different methods to solve 
Freholm integral equations, but none have used coiflets and 
interpolation methods. In this article, we use a scaling 
function interpolation method to solve linear Fredholm 
integral equations of the first kind, and we prove a 
convergence theorem for the solution of Fredholm integral 
equations. We present three examples which have better 
results than previous article. 
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Introduction 

Recently, the study of the Fredholm integral equations 
has attracted the attention of mathematicians and 
researchers in other disciplines such as [6, 7], where 
many problems lead to the solution of integral 
equations [1-3]. Integral equations are also useful in 
many branches of pure mathematics. 

Wavelets have been applied in a wide range of 
engineering and physical disciplines, and are a useful 
tool for scientists and mathematicians. In this paper 
we will find the numerical solution for linear 
Fredholm of the first kind of the form 

b 

f(x) = jk(x,t)y(t)dt, (1) 

a 

where the functions f(x) and k(x,t) are given 
functions over the given interval [a,b] and the square 
respectively, such that the function k(x,t) is called the 
kernel and the unknown function y(t) is to be 
determined. 

Coiflet 

In this section we review briefly on wavelet transform 
and Multiresolution Analysis (MRA) [4]. We first 



define the scaling function (f)(x) and the sequence 
{a p ,p e Z) such that 

</>(x) = Y J Ct^{2 1 x-p) = Y J « Ap (*) ( 2 ) 

p p 

By using this dilation and translation [4], we defined a 
nested sequence spaces [Vj, j e Z} which is called 
MRA of ]} (R) with the following properties 

V^V^JeZ (3) 
V^=f]Vj={0) (4) 

[jVj isdenseinL 2 (i?) (5) 

j 

(/){x) g V, O (/>{2x) e V J+1 (6) 
For the subspace V, is built by 0(2x - p), p e Z then 
V = {(/){x - p),p e Z} And; since VJ, c Vj we have 

p 

In general, 

<j)(x) = X a/(2 j x-p) = ^ a p <j) jp (x) (7) 
p p 

which shows (2) is well defined. Any function 
y(jc)eZ?(7?) can be approximated by scaling 
functions in one of the subspaces in the given nested 
sequence. In fact, for each j we define the orthogonal 
complement subspace W. of V- in the subspace V j+l . 

The orthogonal basis of Wj is denoted by 

y/ jp (x) = y/(2 j x-p), (8) 
and the wavelet function can be obtained by 
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(9) 



Hence the equation (15) becomes: 



More precisely,^, is the subspace generated by } . 

Some interesting properties of scaling and wavelet 
functions make wavelet method more efficient than 
quadrature formula methods, spline approximations 
or other methods in solving Integral equations. A lot 
of computational time and storage capacity can be 
saved since we do not require a huge number of 
arithmetic operations partly due to the following 
properties. 

(i) Vanishing Moments: 

P jcV (x)dx = 0, k = 0, ,m-l. (10) 

J— co 

and in this case the wavelet is said to have a vanishing 
moment of order m, 



(ii) Semiorthogonality: 

Yi. P (*). (*)) = £ ^j, P WVj,k (*) dx = °; ( n ) 

j,k,peZ. 

The set of scaling functions {^„ } are orthogonal at the 
same level n, which means: 

(•CO 

(Kk ( X )A, P ( X )>=L ( X )A, P (x)dx=s kp , 

n, k, peZ 



(12) 



Coiflet (of order L) has more symmetries and is an 
orthogonal multiresolution wavelet system with, 



M 



jx k <f>(x)dx = 0,k=l,2, ,L-l. 



(13) 

lx k y/(x)dx = 0, k = 0,1, L-l. (14) 

Where { M k }are the moments of the scaling functions. 

Scaling Function Interpolation 

The function f(x) can be interpolated by using the 
basis functions in the subspace V ■ as follows. 

f i (x) = Y j a p m j x-p) (15) 

p 

where a p are the coefficients evaluated by using 
equation (12) such that 

a p = (f(x),0 jp (x)} = jf(x)0(2 J x-p)dx. (16) 



/ 1 (*) = Z (f /(^)^( 2; * - P)dx)<f>(2 j x-p). 

p 

On the other hand, one can use sampling values of / at 
certain points to approximate the function /. It has 
been proved in [5], namely, an interpolation theorem 
using coiflet such that if </){x) and y/{x) are 
sufficiently smooth and satisfy the equations (10)-(14) 
and the function /(x)eC t (D) , where Q is a 
bounded open set in R 2 , k > N > 2, j € Z Then, 

f } (x,y) = Wf^^A P m,M 

(x,y)eQ(l7) 
where the index set is 

A = {(p, q)\(supp(0 j p ) ® supp(^ A )) n Q * (/)} ■ 
In addition, the moment M t satisfies 



M,=(cy,/ = 1,2,..., AT-1, 



c = M v 
Then, 



f-f j \\ LHn <c\\f (N) w x (-r.ds) 



1 \N 

'2 j ' 



Where C is a constant depending only on /V, diameter 
of Q and 



(N) 



:- max 

°° (x,y)en,m=Q,...,N 



d N f 



dx m dy 



.(19) 



For the function with one variable, using coiflet (c=0), 
we have 

f i {x) = ^jY j f{^ Up {x),x & [a,bl (20) 

2 p 2 



and 



f-f 1 



L 2 [a,b] 



<c 



co 2 



(21) 
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Where 



:= max 

°° xe(aJ>),m=0,...,N 



dx n 



(x) 



Solutions of Fredholm Integral Equation of 
the First Kind Using Coiflet 

In this section we will apply coiflet and the 
interpolation formula (20) to solve the Fredholm 
integral equation (1). The unknown function y(x) in 
equation (1) can be expressed in term of the scaling 
functions <p j (jc) in the subspace V. such that 



(22) 



Consider the equation (1) and the function y(x) which 
is defined on the interval [a, b] and the scaling 
function <j){x) is supported with (d lf d 2 ) , then we 
have the index: 

A = {2 j b-d 2 ,2 j b-d 2 +l, ,2'a-d^ . This 

coincides with the two dimension index set shown 
above. The reason is as follows. Since the support of (f> 

is (d l , d 2 ) , the nonzero part of 0j (x) is located on 

the interval which satisfies d l <2 J x — p <d 2 . Hence 

d, + p d 7 + p ^ , , , 
1 — <x< — . On the other hand, the 
V V 

unknown function is defined on the interval [a, b]. 

_ . . d, + p , , d 2 + p 

This implies : — < a and b < : — . 

V 2 J V 

Consequently, p<2 J a — d x and 2'b — d 2 < p. This 
proves the index set is shown as above. 

By applying equation (22) into equation (1), we get the 

system, f(x) = j b k(x,t)^a p 0(2 j t - p)dt, (23) 

p 

fix) = ^a p [k(x,ty>(2 j t-p)dt), (24) 

where the coefficients {a p ,pe A] can be evaluated 
by substituting {x p e [a,b], peA} into the system. 

This gives rise to coefficients in (24) and we obtain a 
numerical solution of (1). The error analysis of this 
method is obtained in the following section. 



Error Analysis 

In this section will discuss the convergence rate of our 
method for solving linear Fredholm integral equation 
of the first kind (1) by using coiflets. 

Theoreml. In equation (1), suppose that the function 
k(x,t)eC([a,b] x [a,b]), is positive for all 
x,t e [a,b], Such that 

m l = min{k(x,t), (x,t) e [a, b]x[a, b]}. 



Given f(x) e C[a, b],for j e Z, 

y J (x) = Y J a „0( 2J x ~ P) 



(25) 



is an approximate solution of the equation (1) with the 
coefficients obtained in (24). Then, 

\\e(x)\\ = \\y(x)-y j (x)\\<C(-y (26) 

where, 



u 

\\e{x)\\= \\e{x)\dx 

a 

Proof. Subtracting equation (25) from equation (1) and 
taking the norm for both sides, we get the following 

e(f) = Y. a pK ( f)-y ( f) 

p 

b b 

K = | k(x, t)e(t)dt = J k(x, r)(£ a p (p i p (t) - y{t))dt 

a a P 

b 

II K INI J*(*,f)(X a p( p. p (t) - y(t))dt II (27). 

a P 

b b 

Butfc(x,?) > m x then, j k (x, i)e{t)dt > m^e(t)dt, 



This leads to 



| k(x,t)e(t)dt 



> tnA\e\\ 



Then equation (27) becomes: 
1 



e < 



m. 



£ k(x,t)dt £ (£a/(2 J t- p) - y(t))dt 
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= C, 



where c { = 



[(£a p m J t-p)-y(t))dt 

p 

b 

jk(x,t)dt 



(28) 



By [5], the unknown function y(x) can be interpolated 
by using coiflet such that: 



yJ( x ) = Y J y(^-yp(2 j x-p). 

p 



(29) 



Let t = X in equation (28) then adding and subtracting 

the quantity V y(—)0(2 j t- p) and using the 

p 2 J 

triangle inequality in equation (27), we get the 
following inequality. 



\\e(x)\\ < q 



[(Za p ^(2 J t-p)-y(t) 



f(Zy(^¥( 2 ^-p)-y(0)* 



+ 



(30) 



which equals to the following 



f y(^-)<fi(2't -p)- y(t))dt 



+ 



^y^)-%)\ b J^t-p)dt 



By Theorem 3.1 of [5] and using coiflets, we have 



Iy(Wr-p)-y(0 



1 



^ c o / w (^7)" ( 31 ) 

II I loo O ./ 



Since ^ (y(— ) — a p ) is finite we define it as 
p 2 

2>(^)-a P )=c 2 . 



(32) 



Using the above results and the orthonomality of the 
scaling functions {<p- p (x)} we conclude that 



||e(*)|| < Cl (c \\f m \\(±-) N + c 2 {\V) = c{\y. 



Where c is some constant independent of j. 
Numerical Examples 

In the following examples, we will solve linear 
Fredholm integral equation (1) using coiflet of order 5 
and provide errors between exact solutions and 
numerical solutions at different resolution levels. 
These examples are also presented in [7] by using 
different methods. We will compare our results with 
the results in [7]. 

Example 1 

Consider the integral equation (1) with 
f(x) = e x +(l-e)x-l, 



k(x,t) = 



\t(x-l) t<x 
I jc(/ — 1) X < t 



and the exact solution y(x) = e x . The numerical 
results are shown in Table 1. 

Example 2 

Consider the integral equation 



(1 + x 1 ) 2 -x 



= JV +t 2 ) 2 y(t)dt, 



where the exact solution is y(x) = X and the 
numerical results are shown in Table 2. 

Example 3 

In the following integral equation 



e x+l -\ 



— = \e x 'y{t)dt 
+ 1 J „ 



with the exact solution y(x) = e x . We present the 
numerical results in Table 3. 

In what follows, we present our numerical solutions at 
different levels of resolutions for the above three 
examples. 
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TABLE 1 NUMERICAL RESULTS FOR EXAMPLE 1 



x i 


Exact 




Approximation 


Solution 






solution 


i=-2 


i—1 

1 


i=0 


0.0 


1 


1.00001 


1 


l 


0.1 


1.105171 


1.10516 


1.10517 


1.10517 


0.2 


1.221403 


1.22139 


1.2214 


1.2214 


0.3 


1.349859 


1.34986 


1.34986 


1.34986 


0.4 


1.491825 


1.49189 


1.49183 


1.49183 


0.5 


1.648721 


1.64873 


1.64872 


1.64872 


0.6 


1.822119 


1.82211 


1.82212 


1.82212 


0.7 


2.013753 


2.01377 


2.01375 


2.01375 


0.8 


2.225541 


2.22556 


2.22554 


2.22554 


0.9 


2.459603 


2.45964 


2.45961 


2.4596 


1 


2.718282 


2.7183 


2.71829 


2.71828 


TABLE 2 NUMERICAL RESULTS FOR EXAMPLE 2 


*i 


Exact 
Solution 


Approximation 


Solution 






[ j=- 2 


j=-l 


j=0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


0.1 


0.1 


0.1 


0.1 


0.2 


0.2 


0.2 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


0.3 


0.4 


0.4 


0.4 


0.4 


0.4 


0.5 


0.5 


0.5 


0.5 


0.5 


0.6 


0.6 


0.6 


0.6 


0.6 


0.7 


0.7 


0.7 


0.7 


0.7 


0.8 


0.8 


0.8 


0.8 


0.8 


0.9 


0.9 


0.9 


0.9 


0.9 


1 


1 


1 


1 


1 



TABLE 3 NUMERICAL RESULTS FOR EXAMPLE 3 



Xi 


Exact 
Solution 


Approximation Solution 


j=-2 


j=-l 


j=0 


0.0 


1 


1.00001 


1.00009 


0.99997 


0.1 


1.105171 


1.10516 


1.10517 


1.10517 


0.2 


1.221403 


1.22139 


1.2214 


1.2214 


0.3 


1.349859 


1.34986 


1.34986 


1.34986 


0.4 


1.491825 


1.49189 


1.49183 


1.49183 


0.5 


1.648721 


1.64873 


1.64872 


1.64872 


0.6 


1.822119 


1.82211 


1.82212 


1.82212 


0.7 


2.013753 


2.01377 


2.01375 


2.01375 


0.8 


2.225541 


2.22556 


2.22554 


2.22554 


0.9 


2.459603 


2.45964 


2.45961 


2.4596 


1 


2.718282 


2.7183 


2.71829 


2.71828 



In the following two tables we will compare our 
results with [7] for the first and third examples. The 
results in the second example are close to the perfect 
solution. 

TABLE 4 COMPARISONS OF EXAMPLE 1 WITH OTHER 
METHOD 



The absolute error for example 1 


Our method 


Method in [7] 


*i 


J — 1 


i— 1 
J— 1 


i— n 
J-U 






0.0 


1.0E-5 








4.07E-2 


6E-4 


0.1 


1.1E-5 


1.0E-6 


1.0E-6 


9.99E-2 


1.0E-6 


0.2 


1.3E-5 


3.0E-6 


3.0E-6 


7.7E-5 


3.0E-6 


0.3 


1.0E-6 


1.0E-6 


1.0E-6 


1.41E-1 


1.0E-6 


0.4 


6.5E-5 


5.0E-6 


5.0E-6 


2.5E-5 


5.0E-6 


0.5 


9.0E-6 


1.0E-6 


1.0E-6 


7.71E-4 


1.0E-6 


0.6 


9.0E-6 


1.0E-6 


1.0E-6 


8.4E-4 


1.0E-6 


0.7 


1.7E-5 


3.0E-6 


3.0E-6 


1.52E-3 


3.0E-6 


0.8 


1.9E-5 


1.0E-6 


1.0E-6 


8.71E-4 


1.0E-6 


0.9 


3.7E-5 


7.0E-6 


3.0E-6 


1.60E-2 


7.0E-6 



Where E-n denotes 10 ". 

TABLE 5 COMPARISONS OF EXAMPLE 3 WITH OTHER 
METHOD 



The absolute error for example 3 


Our method 


Method in [7] 


Xi 


j=-2 


j=-l 


j=0 






0.0 


1.0E-5 


9.0E-5 


3.0E-5 


1.11E-3 


5E-5 


0.1 


1.1E-5 


1.0E-6 


1.0E-6 


5.0E-4 


4.1E-5 


0.2 


1.3E-5 


3.0E-6 


3.0E-6 


3.33E-4 


5.7E-5 


0.3 


1.0E-6 


1.0E-6 


1.0E-6 


4.6E-4 


6.9E-5 


0.4 


6.5E-5 


5.0E-6 


5.0E-6 


6.0E-4 


5.5E-4 


0.5 


9.0E-6 


1.0E-6 


1.0E-6 


2.25E-3 


2.4E-4 


0.6 


9.0E-6 


1.0E-6 


1.0E-6 


8.39E-4 


6.9E-5 


0.7 


1.7E-5 


3.0E-6 


3.0E-6 


7.83E-4 


8.7E-5 


0.8 


1.9E-5 


1.0E-6 


1.0E-6 


4.9E-4 


8.8E-4 


0.9 


3.7E-5 


7.0E-6 


3.0E-6 


8.6E-4 


8.7E-5 
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Conclusion 

In this work, we applied scaling function interpolation 
method by using coiflets to solve Fredholm integral 
equations of the first kind. As it is shown in Table 4 
and Table 5, we obtained better results than [7]. Our 
method can be applied to different kinds of integral 
equations, integral-algebraic equations and integral 
equations of two dimensions which will be an efficient 
method for image processing. In addition, we can 
interpolate the given functions in the integral equation. 
This would simplify the calculations in finding 
numerical solutions of integral equations. Using our 
method to solve nonlinear integral equations would be 
interesting as well. 
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